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Abstract 

We consider the large-N Calogero-Marchioro model in two dimensions in the Hamiltonian 
collective field approach based on the 1/N expansion. The Bogomol'nyi limit appears in the 
presence of the harmonic confinement. We investigate density fluctuations around the semiclassical 
uniform solution. The excitation spectrum splits into two branches depending on the value of the 
coupling constant. The ground state exhibits long-range order. 
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I. INTRODUCTION 



The generalized quantum N-body Calogero-Marchioro (CM]model in two and arbitrary 
number of dimensions has been in focus of research recently Q, 0, S, . There are many 
exact solutions and their properties known by now. However, the analytical solution of the 
full many-body CM problem in two dimensions is not tractable. The quantum problem may 
not be solvable even for the ground state |5( for arbitrary values of coupling strenghts of 
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the two and three-body interactions. Thow it is possible to 
radial eigenstates, the set is not complete 

model deserves special attention due to its relation to quantum dots, quantum Ha 
random matrix theory, anyon physics and other systems of physical interest 
Apart of its physical application, we believe that CM model in two dimensions deserves a 
careful field theoretical analysis. In this paper we investigate the large-N CM model in the 
bosonic picture using the large-N collective field technique. 

II. COLLECTIVE FIELD HAMILTONIAN 

The 2D Calogero Hamiltonian describes a system of N non-relativistic particles on a 
plane, interacting via the two-body inverse-square potential and a long-range three-body 
interaction. It is given by 

^ i ijtj \ L t 1 j\ i+j+k I 1 l 3\ I 1 * lfc l ^ i=l 

The particles are confined in a one-body oscillator potential. We choose units such that 
h — m — 1; g and G are in general arbitrary positive dimensionless coupling strenghts of 
the two and the three-body interactions. 

The singularity of the Hamiltonian (1) at = r,,- requires the wave function to have a 
prefactor of Jastrow-type which vanishes for coincident particles. We extract this prefactor 
in the form 

V>(ri, r 2 , r N ) = ]J |r< - r^cj)^, r 2 , r N ) (2) 

i<j 

and obtain the new Hamiltonian 
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acting on the residual wave function <fi. The above Hamiltonian (3) simplifies drastically in 
the limit g = G = X 2 /2 which we use from now on. This is precisely the case when the 
quantum mechanical ground state and an infinite tower of some excited eigenstates may 
be obtained exactly 1]. We are now in a position to introduce the collective field and to 
formulate a non-relativistic quantum field theory which operates in a large N-particle sector 
of the Hilbert space Q|. The collective field is assumed to be the density function 

p(r) = Y,Kr-r l ) , (4) 

i 

where r^'s are positions of N spinless bosonic particles. From the definition (4) it follows 
that the collective field obeys the normalization condition 



J d 2 rp(v) = N . (5) 



Next, we reformulate the differential operators V« in the Hamiltonian (3) in terms of a 
functional differentiation with respect to the collective field p(r) 

= f sk- (6) 

In the large-N limit, the Hamiltonian (3) can be expressed entirely in terms of the collective 
field p(r) and its canonical conjugate 

T(r) = -ij^r ■ (7) 
8p(r) 

The continuum limit of the Hamiltonian (3) can be written as 

H = ^J d 2 rp(r)(Vvr) 2 + y / d 2 rp(v)r 2 + 

+ i\J d 2 r (f J dVp(r)^— ^p(r')^ 7r(r) + i^^- J d 2 r (Ap(r)) 7r(r) . (8) 

In order to re-express all the sums over i and j in the Hamiltonian (3) in terms of the 
collective field p(r), we have been forced to include terms % = j, too. These terms, however, 
are superfluous in the correct form of the Hamiltonian and therefore should be substracted. 
This is the origin of the A dependent part of the last term in the Hamiltonian (8) Q|. The 
Hamiltonian (8) is not hermitian in terms of the collective field p(r) and 7r(r) due to the 



presence of two imaginary terms. In order to obtain the hermitian Hamiltonian, we have to 
rescale the wave functionals using the Jacobian of the transformation from to p(r): 
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The Jacobian J is determined from hermiticity condition and is given by 
A - 2 r ,9 A 



In J = 



J d 2 rp(r)\np(r) + ^J J d 2 rd 2 r'p(r) In |r - i>(r') . 



(9) 



(10) 



The straightforward algebra gives the hermitian Hamiltonian 

2 - A Vp(r 



H= l -jd 2 rp{r)iy^f+ l -Jd 2 rp{v) 



p(r) 
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+ y / d 2 rp(y)v 2 + 
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d 2 rd 2 r'S(r - r')A5(r - r') - XnN6(0) . 



(11) 



The last two singular terms do not give a contribution in the leading order in 1/N expansion 
and should be canceled by the infinite zero-point energy of the collective field p(r). 

III. GROUND STATE AND FLUCTUATIONS 



In order to find the ground state energy of the system, we assume that the corresponding 
collective field is static and has consequently a vanishing momentum 7r(r). Therefore, the 
leading part of the collective Hamiltonian (11) in the 1/N expansion is given by the effective 
potential 

2-A Vp(r) 



Veff = \ I d 2 rp(r) 



P(r) 
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y J d 2 rp(r)i 
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A lower bound on the energy is obtained by performing partial integration in order to rewrite 

V eff as 

r 2 - A Vp(r) 



eff 
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The first term in (13) is positive semidefinite, whereas the rest is constant. The Bogomol'nyi 
bound is saturated by the positive solution po( r ) °f the equation 

2- AVp(r) 



P( r ) 



(14) 



with the ground-state energy equal to 

E = N ( N ~ ^ Xu + Nu . (15) 

This is the exact result [lj . Applying the gradient operator to the integrodifferential equation 
(14), we obtain a differential equation of the Liouville type 

2 — A 

— Alnp(r) - 2vrAp(r) + 2u = . (16) 

It is evident that the character of the solution depends crucially on the value of the parameter 
A. For A < 2, the solution has a Gaussian fall-off at large distances. To our knowledge, its 
analytic form has not been found yet. For special value of A, i.e. A = , the equation (16) 
can be solved exactly by the Gaussian. It is obvious that in the case A / there always 
exist an interesting uniform solution of Eq.(16) given by 

* = (17) 

Needless to say, it's existence is possible on the compact support only. This uniform back- 
ground configuration describes a condensed state of particles. From now on, we will be 
primarly concerned with analysing the spectrum of low-lying excitation around the config- 
uration (17). Performing the 1/N expansion of the collective field p(r) in the form 

p(r) = po + rj(r) , (18) 

where po is the aforementioned uniform ground-state configuration and r/(r) is a small density 
fluctuation around po, we can rewrite the collective Hamiltonian (11) up to quadratic terms 
in 7r(r) and T](r) as 

po f , 2 ,2 , Po f , 2 ( 2 - A V?y(r) f 2 , r - r' V 
H=-jdr{V.) +-Jdr^— -—Xjdr V (r)^-^j + 

+ E + <2 —^- [ [ d 2 rd 2 r'5(r - r')A5(r - r') - \nN5(0) . (19) 



There are no terms linear in rj (r) as we expand around the minimum of the dominant, large- 
N effective potential. Calculation of discarded terms of third order gives an interaction 
between elementary excitations of a Calogero system. Rewriting the second term in (19) we 
obtain 

H = E + ^J rfV(Vvr) 2 + ^£ / d 2 r(V V f+ 



' d 2 rd 2 r'i](r)5(r — r')r/(r') — A 2 7rp / d 2 rd 2 r'r)(r)) In |r — r'\i](r') + 



+ y J d 2 rd 2 r'5{v - r')A<5(r - r) - A7riV5(0) . 

Here, we have employed the identity for the two-dimensional Green's function: 



A In |r — r'l 



r — r 



r — r 



- = 27r5(r - r) . 



(20) 



(21) 



Apart from the long-range two-dimensional Coulomb repulsion (In r) this Hamiltonian con- 
tains the effective hard core, two-dimensional 5-function interaction. As a next step we 
rewrite the density fluctuation , q{v) and conjugate momentum 7r(r) in terms of operators 
a(k) and a^(k) 

^( r ) = ^Jd 2 k (e ikr a(k) + e - ik V(k)) , (22) 



7r r 



J d 2 k (e ikr a(k) - e-^a^k)) , 



2(2tt) 2 ^ 

satisfying bosonic canonical commutation relations 



(k^a^k')! =<J(k-k') 



[a(k),a(k')] = a^(k),a f (k n 



. 



(23) 



(24) 
(25) 



The Hamiltonian becomes 



H = E + J d 2 ke(k)a\k)a(k) + - J d 2 kA(k) (a(k)a(-k) + a\k)a\-k)) + 

+ 5(0) J d 2 ke(k) + J J d 2 rd 2 r'5(v - r')A<5(r - r') - XnN5(0) , (26) 



with 



e(k) 
A(k) 



k 2 /2-A „ . 47r 2 A/V 
+ ——nk + „ 
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k 2 
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'2-A 
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(27) 



(28) 
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It contains non-diagonal term and divergent normal reordering correction to the ground 
state energy, proportional to 5(0), i.e. to the area of the system. To proceed further, we 
have to introduce the new operators 6(k) and tf(k), which are related to a(k) and a^(k) by 



a Bogoliubov transformation 



17] 



a(k) =« k 6(k)+/3 k 6 t (-k) 



(29) 
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a f (k) = ^(k) + /5k&(-k) . (30) 
These new operators satisfy cannonical commutation relations 

[&(k),&t(k')]=5(k-k') , (31) 

[6(k),6(k')]=[6 t (k),6t( k ')]=0 . (32) 

provided that 

a£-# = l. (33) 

The coefficients «k and {3^ of the Bogoliubov transformation are determined by requirement 

that non-diagonal term of the Hamiltonian (26) expressed in terms of operators &(k) and 
&t(k) vanish and this condition leads to 

2 e(k)+o;(k) 

" k = Mk) ' (34) 

2 e(k) - u(k) 

^ = Mk) • (35) 

Then the Hamiltonian is diagonalized and reads 

w(k) , 



H = E + J d 2 kou(k)tf(k)b(k) + 5(0) | d 2 A; 

+ y J d 2 rd 2 r'5(r - r')A5(r - r') - XnN5(0) , (36) 

where u(k) represents the energy spectrum of the physical fluctuations 

2-A 



u(k) = ^(k) - A*(k 
The ground state |0) is defined by 



-k z + 2u 



(37) 



6(k)|0) = (38) 
and &t(k) is the creation operator of the low- lying excitation with momentum k. 



IV. DISPERSION LAW AND DIVERGENCIES 

We note that there are two qualitatively different dispersion laws, depending on the value 
of the coupling A. For A < 2, the function o>(k) has one stationary point at k = (absolute 
minimume) 

uj(k) = ^-^>k 2 + 2u, A < 2 . (39) 
7 



For A > 2, the function u;(k) has one stationary point at k = and infinite number of zero 
points located at the circle with critical radius k c = \fj 



_Suj_. 
A-2" 



W (k) 



^k 2 + 2lo , |k| < k c 
^k 2 - 2uj , |k| > k c 



, A > 2 . (40) 
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Thus, there is a discontinuity in slope at k — k c . The minimum of oo(k) touches zero 
at |k| = k c . The appearance of such gapless excitation with nonzero momentum can be 
considered as an indication of instability of the system to phase separation, or collaps. This 
result clearly shows the breakdown of the Bogoliubov approximation near the critical radius 
|k| = k c . The coefficients « k and /3 k diverge there, indicating that collective-field dispersion 
cannot be trusted for large values of the quasiparticle momentum. We note that in both 
phases there is the plasmon gap at k = 0, generated by the induced two dimensional long- 
range Coulombic interactions. As long as A > 2, the low- momentum excitation spectrum 
(40) produces a minimum between the 'plasmon' frequency 2u and the 'free-particle' limit 
lu ~ k 2 /2m*, where m* is the effective mass of the quasi-particle: 

2 

m* = — —m . (41) 

I — I 

Notice that the quasiparticle effective mass m* increases with A exhibiting a divergence at 
a critical value A = 2, after which it decreases with increasing A. The above dispersion 
relations have interesting interpretation in terms of the effective, Fermi-sea excitations. A 
single hole-like excitation is obtained by moving a state from the Fermi sea (|k| < k c ) to 
the Fermi surface (|k| = /c c ), while a single-particle-like excitation is obtained by moving a 
state from the Fermi surface to a level above (|k| > k c ). There is nothing strange about 
the presence of such fermionic features in our genuine bosonic picture. After all, we have 
assumed that our wavefunctions incorporate hard-core condition in the form of the Jastrow 
prefactor. 

Now we consider the problem of singular contributions present in the Hamiltonian (36). 
Using the integral representation of the Dirac delta function we can rewrite the singular part 
of the Hamiltonian (36) as: 

2- A 



H 



sing. — g 



J d 2 rA5{r) 



r=0 



- XnN5(r) 



r=0 



A-2 A r „ XttN 



8 ^/^-(S/^- (42) 



where A denotes the infinite area of the system 

A= [ d 2 r = (2vr) 2 5(k)| k=0 . (43) 



Since the total number of particles N and the area A are interrelated by relation (17) 



N uj 

A ATT 



the singular part (39) reduces to 



J^sing. — n 



d 2 k(^-^k 2 -2oj) . (44) 
k=o - ' V 4 / 

In the case of the first phase A < 2, it is evident that the singular contribution H S i ng . 
completely cancels the divergent normal reordering correction. There is no correction to the 
ground state energy due to quadratic fluctuations. Naively, in the phase A > 2 there is no 
complete cancelation and the correction to the ground state energy is found to be: 

A- 2 r°° 



A£o = 5(k)| k=0 '—^ I d 2 k(k 2 - k 2 c ) . (45) 

4 J\k\=k c 

However, taking into account the relation (17) and the definition of k c we see that in our 
approach the sector |k| > k c is absent because of large-N limit, i.e. k 2 ~ po ~~ > oo. There is 
no correction to the ground state energy again. The finite Hamiltonian of quantum collective 
excitations reduces to 

H = 1 -^-—^- \oj + Nu + J d 2 k^^k 2 b\k)b(k) + 2ujjV , (46) 

where the number operator j\f labels the quanta of oscillator excitations: 

N = J d 2 kb\k)b{k) . (47) 

The excitation spectrum of density fluctuations around the semiclassical solution pn shows 
roughly the infinite tower structure of eigenvalues, separated by the spacing 2uj Each 
level is smeared by small, continuous contributions from kinetic energies of particles (A < 2) 
or holes (A > 2). 

V. WAVEFUNCTIONALS 

Now we are in a position to find the collective-field vacuum wave functional. With the 
definition (38) and the Bogoliubov transformation (29) and (30), we get 

6(k)|0> = (a k a(k) - /W(-k))|0) = . (48) 



Rewriting the bosonic operators a(k) in terms of the density fluctuation rj(r) and conjugate 
momentum 7r(r) the vacuum condition (48) becomes 

6 



I 



d 2 re i]ir 



(« k - fa)v(r) + 2(27r) 2 po(ttk + fa) 



5i](r) 



|0> = . 



(49) 



It is evident that solution of this linear functional derivative equation is given by the Gaussian 
ansatz: 

|0) = el d2rd2r ' r l(r)K(r-r')ri(r') 



with translationally invariant kernel K(r). Using Fourier transform 

K(r - r') = J d 2 ke lk{r ~ r,) K(k) 

and combining relations (49), (50) and (51), we get 

J d 2 re lkr V (r) [(a k - fa) + 4p (2vr) 4 (a k + fa)K(k) 

Since this equation holds for any k, the kernel K(k) is given by 

1 a k - fa 



(51) 



(52) 



K(k) 



1 



A - 2 2u 
~~ 4 



(53) 



4(2vr) 4 po« k + /3 k 2(2tt) 2 P o 

This result enables us to reconstruct the Schrodinger wave function ■0(r 1 ,r 2 , ■■■,r N ) for the 
ground-state of the N-particle system. The ground state functional is given by the 



*[p] = nir,-i;|\H| 



(54) 



Here, the prefactor is present owing to the extraction (2). The Jacobian of the transformation 
from ri to p(r) rescales wavefunctional by the factor J^. Expanding the Jacobian to the 
quadratic term in r\ and using relations (10), (17) and the Bogomol'ny equation for p (14), 
we are left with 



i>j 



~ J d 2 rr 2 r](r) - ^ J J d 2 rd 2 r'r)(r) In |r - r'\rj(r')+ 



+ 



2 - A 



— J d 2 rrf{r) + J J d 2 rd 2 r' V {r)K{r - r') V (r') 



8p, 

Next, we employ the Fourier representation of the In kernel: 

gik(r-r') 



In |r 



d z k- 



k 2 



(55) 



(56) 
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so that the vacuum functional \I/ can be rewritten in the A < 2 phase as 

* = n|r i -r i |VS/« ft -'V'> . 



(57) 



If we substitute equation (4) in (57), we obtain the well known ground-state wavefunction 
for N-particles j]]: 

*(r 1 ,r 2 ,...,r JV ) = nir<-r i | A e-^5:i r ? . (58) 



In order to find the vacuum wave-functional in the A > 2 phase, we rearrange the contribu- 
tions from integrations in (55) and obtain the following expression 



i>j 



r r x 



x exp 



— j [ d 2 rd 2 r'r](v)ri(r') ( [ d 2 k 

2,71 J J \J\k>kc 



gik(r-r') 



gik(r-r') ~ 



d 2 k , 

|k|>fc c 



(59) 



'|k>/r |k| 2 

We get the same vacuum wave-functional as in the A < 2 phase (57), since k c — > oo in the 
large-N limit. The unrenormalized wave function of the first excited state may be represented 



by 



^(k) = 6 f (k)|0) . 



(60) 



Using the complex conjugate of relation (48) it is easily shown that the phonon state is 
mearly the coordinate multiplied by the Gaussian form: 

2 



V(k) = —fa-fo) I d 2 n7(r)e ikr |0> . 



(61) 

When expressed in configuration space coordinates, these excited states are just the functions 
used by Feynman to describe phonon and roton states [l^ : 



V(k)~£e* kr <|0> . 



(62) 



We emphasize here that the one-particle excited states have vanishing norm at |k| = k c . 
This is the consequence of the fact that |k| = k c is the singular point of the Bogoliubov 
transformation (34) and (35) (a(k c ) = /3(k c ) = oo). However, these large momentum 
excitations are beyond the scope of our approximation. 
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VI. CORRELATION FUNCTION AND STATIC STRUCTURE FACTOR 



Next, we study the density- density correlation function and the static structure fac- 
tor in the 2D CM model. The analysis is similar to the one dimensional case (2^]. The 
density— density correlation function G(r) is defined by 

(0|p(r)p(r')|0> = pi + (0\r)(r)rj(r f )\0) = & + G(r - r') . (63) 

The relations (22), (29) and (30) enable us to evaluate the function G(r) 

G( T ) = -A- [ dV^j- . (64) 
V ; (2tt) 2 J 2u(k) v ; 

The static structure factor S(k) is defined as the Fourier transform of G(r): 

This is the familiar Feynman-Bijl relation. By expanding the static structure factor S(\t) in 
powers of k 2 /p up to the cubic terms, we easily obtain 

S(k) k2 (l T ( 2 ~ X \^ 1 ( 2 - A ) 2 l^ T ^ (66) 
Mkj " 4Att y- T SAvrpo k + (SAvrpo) 2 T "') { } 

where the upper and lower sign refers to the A < 2 and A > 2 phase, respectively. We are 
now in position to compare the collective field theory structure factor with the exact one, 
known from the random matrix theory for special value A = 1 0, 0> 13 Namely, it is 
easily seen, that the squared value of the ground-state wavefunction (58) is proportional to 
the joint probability density function of the eigenvalues of complex matrices, provided one 
sets A equal to 1. By expanding the exact static structure factor S\=i(k): 



S x=1 (k) = Po y-e~^j , (67) 

we obtain 

S x=1 (k)=p [- — - + — . (68) 

\4:Trp 2(ATrp y 6(4vrp ) d J 

Comparing relations (66) and (68), one immediately sees that they agree up to forth order 
for |k| close to 0. This indicates that our method is meaningful as long as we restrict 
ourselves to the wavelengths of density fluctuations much longer than the average inter- 
particle separation l/^/po. Having found the static structure factor, we now turn to the 
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calculation of the correlation function G(r). Using the definition (64) and expansion (66) 
we find that 

The asymptotic form of G(r) at large values of |r| gets only a contribution from the small 
|k| region in the integral (69). Using the definition of the Bessel functions of the first kind 
of order zero, it can be seen that the correlation function shows long-range order at large 
distances for any value of A ^ 1: 

GW ~ "8A^ 1 1 <** J ° [X) ± 1 1 J " (X) + - (7 °» 

VII. CONCLUSION 

In this paper we have formulated a quantum field theory of two-dimensional CM model 
in terms of collective-field variables in the large-N limit. In this semiclassical approach, we 
have determined the ground-state energy and the corresponding uniform density solution. 
Using the method of the Bogoliubov approximation, we have found the spectrum of low-lying 
excitations around the ground state for any value of the coupling constant A. We have shown 
that the spectrum splits into two branches depending on the sign of A — 2. We have proposed 
the interpretation of the spectrum in terms of quasi-particle (A < 2) and quasi-hole (A > 2) 
excitations in the effective Fermi sea, defined by the Fermi momentum k c . We have also 
shown that the ground-state energy is not affected by the next to leading order corrections 
steaming from the quantum collective-field fluctuations. Furthermore, we have obtained the 
density-density correlation function and the static structure factor which shows long-range 
order. Finally, we believe that the exact u;(k) does not vanish at nonzero value of |k| for 
A > 2. The correct dispersion probably has a roton-like minimum at |k| = k c . In order to 
test the existence of roton dip in the spectrum of the excitation one has to go beyond the 
quadratic approximation. We hope to report on this issues in a separate publication. 



-k 2 T 



(2 -A) 



4Att 32A% 2 po 



k 4 + ... 



(69) 
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